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Voyager images of the satellites of Saturn and Uranus have shown that these bodies are charac- 
terized by remarkable diversity and surprisingly complex geologic histories. Despite their small 
sizes, a number of the satellites show unambiguous evidence for resurfacing. Satellites that 
have clearly undergone at least one episode of resurfacing include Enceladus, Dione, Miranda, 
Ariel, and Titania [1,2]. Less compelling evidence has also been presented for resurfacing on 
Mimas and Rhea. In the case of Enceladus and perhaps Ariel, tidal dissipation may have been 
largely responsible for heating the satellites. The energy source that apparently led to melt- 
ing in the other satellites, however, is less clear. For all the satellites, long-lived radiogenic 
heating is entirely inadequate to produce melting. We have therefore turned to investigation 
of accretional heating as a possible energy source. Numerous models of accretional heating of 
planetary bodies have been developed previously, especially for the terrestrial planets. All such 
models suffer from a need to select preferred values of critical model parameters in order to 
obtain meaningful results. Our goal here has been to develop a detailed model for the heating of 
these small satellites, and then to explore the consequences of variations in the free parameters 
in the model. Specifically, we attempt to determine for what range of conditions melting will 
occur in these satellites. Along with varying a number of model parameters, we also consider 
the important effects of inclusion of small amounts of ammonia and (in the case of the uranian 
satellites) methane in the system. 

We begin with a swarm of bodies from which the satellite will accrete. We allow accretion to 
take place gravitationally, calculating the amount of heat deposited, the depth over which it 
is deposited, the subsequent transport of heat within the satellite, the rate of growth of the 
satellite, and the satellite’s evolving internal thermal profile. We consider two reservoirs of 
particles from which a satellite may accrete: bodies in orbit about the planet (cis-planetary 
debris) and bodies in orbit about the sun (trans-planetary debris). The initial surface density 
a c of the cis-planetary particle disk is obtained by distributing the present mass of the satellite 
over an annular region centered on the planet, with inner and outer radii at the distances at 
which the gravitational attraction of the present satellite is balanced by the attractions of its 
inner and outer nearest neighbors, respectively. The surface density at of the trans-planetary 
disk is determined by similarly distributing a small fraction /3 of the planet’s mass over a sun- 
centered annular region. The impacting particles are assumed to have a density equal to that 
of the present satellite, and to have a size distribution of the form: 
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where N(r p ) is the number flux of particles of radius r py and a is a constant. 
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In calculating the accretion rate and the impact velocity of the impactors from both particle 
disks, we follow the formulation of Safronov [3]. We use the Safronov parameter 9 C for the 
cis-planetary disk, defined by 

6 C = GM t /R t u c 2 (2) 

where G is the gravitational constant, M $ is the mass of the satellite, R s is the radius of the 
satellite, and u c is the velocity of an impactor relative to the circular orbit. The Safronov 
parameter is related to the eccentricity of the particles in the disk, and is a free parameter in 
the model. It then follows that the impact velocity v c is given by 


v c — \J + 2GM & /R b 


(3) 
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and that the mass flux M e is given by 
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where P s is the satellite’s orbital period. Following Kaula [4], we also define similar expressions 
for vt and Mt, the impact velocity and mass flux for debris from the trans-planet ary disk. 

Having developed expressions for the size distribution, impact velocity, and mass flux of particles 
from both sources, we then calculate deposition of thermal energy by the impact process. We 
assume that some fraction f) of an impactor’s kinetic energy is deposited beneath the satellite’s 
surface; r) is another free parameter in the model. When an impact takes place, a roughly 
hemispheric shock wave centered on the impact point is generated. As it passes through the 
material of the satellite, energy is deposited in that material. The specific energy deposited 
is largest near the impact, and decreases monotonically with distance from the impact point. 
We have solved the Rankine-Hugoniot equations for the deposition of specific energy by shock 
passage, using published experimental data on shock propagation through H2O ice. The specific 
energy e deposited as a function of distance r from the impact point is taken to be 

e = e 0 (r i /r)' 1 (5) 

where r s - is the impactor radius. Solution of the Rankine-Hugoniot equations for ice show that 
an expression of this form is appropriate, and give a value of 7 = 3.1. Energy deposition is 
described by eq. (5) out to the point at which the shock pressure drops below the Hugoniot 
elastic limit for ice [5], and beyond that point is assumed to be zero. Impacts are assumed to 
be uniformly distributed over the satellite. Energy deposition is integrated over depth and over 
the particle size distribution for both impactor sources, so that the total energy deposition as 
a function of depth is obtained. 

The satellite is allowed to grow, increasing in radius as it depletes the particle disks from 
which it is growing. Simultaneously, we solve the heat conduction equation for the growing 
satellite, allowing heat to be transported in response to internal thermal gradients. The thermal 
parameters used are those appropriate for H2O ice. Near the surface, heat transport will also 
be significantly influenced by the physical mixing that takes place due to impacts, introducing 
an additional effective thermal diffusivity. The effective diffusivity due to impact mixing at a 
depth z is given by 

K{z) - In*z 2 (6) 

where h z is the number of times a layer of depth z is turned over by impacts per unit time. We 
have developed an expression for K(z) based on the regolith mixing model of Gault et a 1. [6]. 
Growth of the satellite takes place until the source of impactors is exhausted. The free param- 
eters in the model are then 0 (the fraction of the planet’s mass in the trans-planetary disk) , 
9 C and $i (Safronov parameters for impactors in both disks) , a (the particle size distribution 
exponent), and fj (the impactor kinetic energy partitioning coefficient). 

We have performed two sets of calculations. In the first, we have done a large number of 
calculations of accretional heating of Saturn’s satellite Rhea, allowing for large variations in all 
the paramters of interest. A sample result is given in Figure 1, for 0 = 0.001, 6 e = 4, 9% = 4, 
a = 3.0, fj = 0.2, and 7 = 3.1. Reasonable variations in 0 have a negligible effect on the results. 
Even placing 10% of the mass of Saturn in the trans-planetary nebula has little effect; the reason 
is that the timescale for accretion from the cis-planetary disk is very short in comparison to 
that for the trans-planetary disk. Variations in 9 C are important. For large values of 9 C (i.e., 
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for low impactor eccentricities) , accretion times are extremely short; for example, ~ 200 yr for 
$ e = 50. Impact velocities are also low, so energy is deposited only very near the surface. In 
contrast, low values of 9 e le&d to longer accretion times 5000 yr for 6 C = 2), and deeper 
deposition of energy. Variations in a are also important; smaller values of a mean that more 
of the mass flux is in large impact ors and increases the amount of energy deposited at depth. 
As one would expect, variations in rj have a significant effect on the peak internal temperature 
achieved. 

In the second set of calculations, we have calculated accretional temperature profiles for the 
inner major satellites of Saturn (Mimas - Rhea) and for all the major satellites of Uranus. The 
maximum temperatures achieved are: Mimas 83 K, Enceladus 93 K, Tethys 145 K, Dione 158 
K, Rhea 195 K, Miranda 82 K, Ariel 184 K, Umbriel 160 K, Titania 214 K, Oberon 199 K. 
Given the uncertainties in the input parameters, these temperatures may be uncertain by as 
much as a few tens of degrees. It is clear that temperatures will never exceed 273 K, so melting 
would not be expected if the satellites were made of pure H 2 O ice. However, it is probable that 
the saturnian satellites contain some NH3, and it is possible that the uranian satellites also 
contain some CH 4 . A peritectic melt will form in the H2O - NH3 system at a temperature of 
173 K, and it has been suggested that ice containing CH4 could be substantially mobilized at 
temperatures as low as 100 K [7]. Given the uncertainty in the calculations and the likelihood 
of NHs in the saturnian satellites and CH4 in the uranian satellites, these calculations suggest 
that, among these satellites, only Mimas and Enceladus cannot have undergone resurfacing as 
a result of accretional heating. 
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Figure 1 — Calculated thermal profiles within Rhea during accretion. Curves are plotted at 
equal satellite mass intervals. Total accretion time is ~ 2000 yr. 
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